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Abstract 

A new simple method to measure the spatial distribution of the electric field 
in the plasma sheath is proposed. The method is based on the experimental 
investigation of vertical oscillations of a single particle in the sheath of a low- 
pressure radio-frequency discharge. It is shown that the oscillations become 
strongly nonlinear and secondary harmonics are generated as the amplitude 
increases. The theory of anharmonic oscillations provides a good qualitative 
description of the data and gives estimates for the first two anharmonic terms 

in an expansion of the sheath potential around the particle equilibrium. 
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There are various experimental methods on the plasma sheath diagnostic in low-pressure 
gas discharges (see, e.g., and references therein). Most of the methods are based on 
optical emission spectroscopy and temporally resolved probe measurements and are rather 
complicated technically. We propose an "alternative" simple method to measure the spa- 
tial distribution of the electric field in the plasma sheath based on investigation of "large- 
amplitude" vertical oscillations of micron sized particles. 

So far, investigation of the particle oscillations in the sheath of a radio-frequency (rf) 
discharge 0-0] was one of the methods used to determine the particle charge |]8HlO|. This 
method is based on the assumption that the vertical distribution of the sheath potential 
is known and can be approximated by a parabolic profile |Tl|, so that the oscillations are 



always harmonic. This assumption is reasonable for sufficiently high gas pressure, when the 



mean free path of ions is much less than the vertical spatial scale of the sheath field |jTT| : For 
example, the model of a parabolic rf sheath in an argon plasma can be applied for pressures 
above ~ 20 Pa. But for p ^ 10 Pa the sheath profile might deviate strongly from a parabolic 
one. In this case, the oscillations are harmonic only if their amplitude is very small. 

In this Letter, we report experimental and theoretical investigations of nonlinear (large- 
amplitude) vertical oscillations of a single particle in the sheath of a low-pressure rf discharge. 
We compare the experimental results with the theory of anharmonic oscillations and estimate 
the first two anharmonic terms in an expansion of the sheath potential around the particle 
equilibrium height. Thus, the study of large-amplitude particle oscillations allows us to 
reconstruct the distribution of the vertical electric field in wide region of the sheath. 

Experimental setup and results — The experiments are performed in a so-called Gaseous 
Electronic Conference (GEC) rf reference cell [|1^ with the electrode system modified as 
shown in Fig. |I]. The lower aluminum electrode, 230 mm in diameter, is capacitively coupled 
to a rf generator providing a peak-to-peak voltage ~ 70 V (rf power ^ 2 W). A ring-shaped 
upper grounded electrode with an outer diameter of 200 mm and a thickness of 10 mm is 
located 50 mm above the lower electrode. A 2 mm thick copper ring with a 40 mm inner 
diameter is placed on the lower electrode to confine the particle above the center. Argon 



gas at a pressure 0.5 Pa is used for the discharge. For this pressure, a self-bias voltage 
~ —33 V on the lower electrode provides the sheath with a (visible) thickness of ^ 20 mm 
and causes a spherical polystyrene particle (density 1.05 g/cm'^, diameter 7.6 ± 0.1 fim) to 
levitate ~ 14 mm above the electrode. 

At a height 8 mm above the lower electrode a horizontal wire of 0.5 mm diameter and 
80 mm length is placed directly below the levitated particle. A low-frequency sinusoidal 
voltage of amplitude Uex is applied to this wire to excite the vertical oscillations of the 
particle. The frequency uj/2tt of the voltage could be varied from 0.1 Hz up to 40 Hz and 
back in 0.1 Hz steps. To determine the oscillation amplitude at each step, we illuminated 
the particle with a vertical laser sheet of ~ 140 fim thickness and imaged it from the side 
by an external digital CCD camera (max. field rate is 50 fields/s, resolution is 768 x 576 
pixels) . 

For relatively small values of the excitation voltage {Ucx ~ 30 mV), we observed the usual 
linear (harmonic) oscillations: The dependence of the oscillation amplitude on the frequency 
has a single symmetric narrow maximum at the resonance frequency u;o/27r ^ 17.0 Hz (pri- 
mary resonance). However, for larger f/gx the dependence clearly exhibits nonlinearity of 
the oscillations (see, e.g., |T^). Figure |^ shows that the maximum of the primary reso- 
nance shifts to the left, as the excitation voltage increases. When Ue^ exceeds some critical 
value, a hysteresis appears. To illustrate this phenomenon, we can compare Figs. |^(a) and 
11(b) (f/cx = 100 mV and 200 mV). When the frequency increases [Fig. |(a)], the amplitude 
grows continuously until a certain frequency {u!+) is reached, then it "jumps" upward and 
thereafter decreases monotonically as u is further increased. On the "way back" [Fig. 0(b)], 
the amplitude grows continuously until another frequency (uj-) is reached. Then it sponta- 
neously "jumps" downward and continuous to decrease slowly as the frequency is reduced 
further. We found that is always larger, than uj_. Outside this "hysteresis zone", the 
experimental points obtained for increasing and decreasing frequencies practically overlap. 
The shift of the maximum, as well as the width of the "hysteresis zone" rapidly increase 
with t/cx- 
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Another peculiarity of the oscillations are the secondary resonances - two maxima at 
uj ^ |u;o and uj ^ 2ujq, which are observed for sufficiently high excitation voltage. Figure |^ 
shows the peak in the vicinity uj ~ \uJo (superharmonic resonance) for different values of [/ex- 
We see that the dependence has a finite offset with some slope (marked by the dashed line), 
which represents the low-frequency tail of the primary resonance. The shape of the peaks 
resembles that of the primary resonance: As [/ex increases, the maximum shifts to lower uj 
and the hysteresis appears. In contrast, for the subharmonic resonance at ~ 2cJo (Fig. g) 
the dependence is qualitatively different. When f/ex is below a certain threshold, we see only 
the offset - the high-frequency tail of the primary resonance - but when f/ex exceeds this 
threshold, the amplitude increases abruptly. When oj is increased [Fig. H(a)], the amplitude 
decreases very slowly along the offset line up to some point. Then the amplitude jumps 
upward and subsequently decreases rapidly until the resonance branch crosses the offset 
line. At higher frequencies the dependence again follows the offset. When oj is decreased 
[Fig. g(b)], the amplitude varies approximately as in the case of increasing frequency, but 
there is no readjustment downward anywhere. The amplitude increases until the particle is 
ejected from the sheath. Both the height and the width of the peak rapidly increase with 

t/ex. 

It is worth noting that the amplitude of oscillations both for primary and secondary 
resonances is always much smaller than the sheath thickness (~ 20 mm). 

Theory and discussion — The particle oscillations are determined by the spatial distribu- 
tion of the electrostatic potential, ^(z), in the sheath. In the absence of excitation, a particle 
is levitated in the minimum of the potential well (at the height 2; = 0), where the gravita- 
tional force, M(7, is balanced by the electrostatic force, Q(0)£'(0) = QqEq. We assume that 
the equilibrium particle charge, Q{z)^ changes weakly with z [compared to the electric field 
of the sheath, E{zy\, so that we can set Q ^ Qo = const. The particle collision experiments 



Ijl in Ar discharge at pressures p ~ 1 Pa confirm this assumption: Typically, the spatial 
scale of the field change is approximately one order smaller than that of the charge change 
i TSj . The electrostatic energy of the particle, lA = Q(f), can be expanded around 2; = in the 



U{z) = U'oz + ^U'^z-" + \u;'z^ + ^W^^' + 0{z^). (1) 

Using the condition of equilibrium Mg = QEq = —Uq and the definition of the resonance 
frequency, McUq = —QEq =Uq, we rewrite it as follows: 

U{z) ~ M (-gz + \ulz^ + ^ i^^4^ ^ (2) 

where a = Uq /2M and (3 = Uq^^ /6M are the anharmonic coefficients. We keep only the 
first two anharmonic coefficients, because in practice they are sufficient to describe all major 



peculiarities of a nonlinear oscillator |13,16 



Introducing the damping rate due to neutral gas friction, 7, and using Eq. (Q), we can 
write the following equation for the particle oscillation: 

F(z) 

z + 2'jz + lUqZ + az^ + Pz^ = cos cut. (3) 

Here F is the excitation force acting from the wire below the particle. This force presumably 
depends on the vertical position of the particle, and we expand it up to quadratic terms, 

F(^)~Fo + FoZ + iFoV_ (4) 

First, let us consider the primary resonance u; ~ cjq- Nonlinear effects are stronger 
when the oscillation amplitude is larger, and therefore we are interested in a narrow region 
e = uj — ujq around the resonance frequency, where |e| <^ ujq. This approach is normally 
valid in the limit 7 <^ cuq and is justified for our case: From the expression for the neutral 



friction force [|T^ we evaluate 7/27r ~ 0.068 Hz which is much less than the measured 



tUo/Svr 17.0 Hz. Then, we obtain from Eq. (|^), 

z + 27i; + cUqZ + az'^ + jSz^ = tUgA cos(co'o + e)t, (5) 

where A = Fq/MuJq = kUcx and k is a "scale factor" of the oscillation amplitude, normalized 
to the excitation voltage. In Eq. (1) we omit the terms of higher order in F{z), because for 



the primary resonance they are not important when the inequahty A|Fq/Fo| <^ 1 is satisfied 
(this condition always holds in our case). In accordance with the method of successive 
approximations for anharmonic oscillations 0,|l^], the first term in the solution of Eq. 
is z{t) ~ a cos (cuq + e) t, and the dependence of the amplitude a on e and A (frequency- 
response equation) is given by 

aM(e-«:a^)^ + 7^] = H'A^ (6) 

where the constant k characterizes a nonlinear shift of the primary resonance frequency, 

^ 3^ _ ^ 

The least-squares fit of the experimental points for the low-amplitude (linear) primary reso- 
nance gives us the following values: the resonance frequency ujq/2tx = 17.0 Hz, the damping 
rate 7/2?? = 0.067 Hz, and the amplitude scale factor k = 0.042 mm/V (note that the fitted 
damping rate is very close to 7/2?? = 0.068 Hz evaluated from the theory [0). Using these 
values, we fit the points for the nonlinear oscillations in Fig. ^ to the function a(e. A) from 
Eq. (^. The solid lines represent the least-squares fit with K/27r = —0.96 Hz/mm^. If each 
curve is fitted independently, this value varies within 5%. The hysteresis appears on the 
curves, when A exceeds the critical value: A^^ = ■^^{'-f^ /\k,\uJq) ~ (2.4 x 10~^ mm)^, which 
corresponds to f/ex — 60 mV. 

For the superharmonic resonance cu ~ ^ujq, we set e = — ^Uq. The solution of Eq. (^ 
is approximately a sum of the first- and second-order terms |T^, z{t) ^ zi{t) + Z2{t). The 
first (linear) term is Zi(t) = ai cos(|co'o + e)t, where 

ai ~ I (1 + |e/cuo) A, (8) 

represents the amplitude offset (dashed line in Fig. ||), and Z2 describes the resonance peak 
due to nonlinearity. Substituting the expression for zi{t) together with Eq. (|) in Eq. (^, 
we get the resonance terms zf oc cos(ciJo + 2e)t in the resulting equation (see Ref. [0). 
Retaining these terms, we obtain. 



Z2 + 27i2 + ujIz2 + azl + I3zl 



9«(^/2) 



COS (cuo + 2e) t, 



(9) 



where a 



(^/2) 



a 



C(Jq/£ and £ = Fq/Fq < is the spatial scale of change of the excitation 
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force. In Eq. (P) we omit the small terms 0{A^/l^). Formally, the derived equation for 
the superharmonic resonance is similar to that for the primary resonance [see Eq. (^], and 
therefore shapes of the corresponding curves are qualitatively the same. Assuming Z2{t) ~ 
a2 cos (tuo + 2e) t, we obtain the following frequency-response equation for the superharmonic 



resonance: 



(2e - /tag) +7 
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Hu/2) 



(10) 



^1 V ujQ 

The independent least-squares fit of the experimental points in Fig. |^ using Eq. (|10D gives 
for each curve values of ujq and 7, which are nearly the same as those for the primary 
resonance (deviation within ~ 0.5%), K/27r = — 1.0 ± 5% Hz/mm^, and |a;(a;/2)|/(27r)^ = 
32 ± 10% Hz^/mm. The critical value of A for the hysteresis: A'^^ = ^(7^c<;o/|fi;|«^(^/2)) — 
(4.5 X 10^^ mm)^, which corresponds to Ucx — 3.5 V. 

Next, for the subharmonic resonance u ~ 2ujq, we set e = u; — 2uJo. The first-order term 
is zi{t) = —ai cos(2co'o + e)^ with the amplitude offset 

ai ~ i (1 - le/uo) A. (11) 

Retaining the necessary terms (see Ref. |ITB|), we derive from Eq. (|^) the equation for Z2, 

(12) 



Z2 + 272:2 + oJq 1 - I (ai^2Lo)A/ujl^ cos {2ujq + e)t Z2 + az\ + = 



where a(2uj) = a + fi^o/^- \h.dX Eq. (|T2[) is an equation of a nonlinear parametric 

oscillator, and therefore the subharmonic resonance increases due to a parametric instability 
H. Putting Z2{t) ~ 02 cos(u;o + \^ti we get the frequency-response equation, 
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[ie - fi:a2)2 + 7^ 
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A'al 



(13) 



Solution of Eq. (|I^) has the following peculiarities |]T^ : If A is less than the threshold value 
At\i = 67u;o/|a(2(^)|, then Eq. (|13D has only the zero solution a2 = 0, and thus just the 



offset ai{uj) can be observed in experiments. For sufficiently high Ucx {A > Ath), a non-zero 



solution exists for ebi < e < eb2, where ebi,2 = T'^l y A'^ / — 1. If the frequency increases, 
02 jumps upward from zero at ebi and then continuously decreases to zero at eb2- On the 
"way back" , 02 starts from zero at eb2 and then grows monotonically. The least-squares fit 
of the points in Fig. ^ (f/ex = 2.0 V) with Eq. ([T3|) gives us k/2'k = —0.95 Hz/mm^ and 
|c^(2a;) |/(27r)^ = 320 Hz^/mm. Thus, — 0.02 mm, which corresponds to Ucx — 0.5 V. The 
fit for the center e = yields 34.85 Hz, which is ~ 2% larger than the double value of the 
resonance frequency obtained for the primary and superharmonic resonance. 

Now we can finally reconstruct the sheath potential, (piz). Using the fitted values 
|a(<^/2)| = \a-lujl/£\ - (27r)2x30 HzVmm and |a(2^)| = + ^ (27r)2x320 HzVmm, 

we evaluate the first anharmonic coefficient and the spatial scale of the excitation force: 
a/(27r)^ ~ —130 Hz^/mm and £ ^ —2 mm. Substituting this value for a and putting 
/t/27r = —0.95 Hz/mm in Eq. we estimate the second anharmonic coefficient: 

(3/{2'kY — 20 Hz^/mm^. Then, using the series Eq. (Q) we evaluate the potential energy of 
a particle in the sheath field {z in mm). 



We see that in our case at least the first anharmonic term results in a significant deviation 
of the sheath potential from a parabolic one, and therefore the oscillations become strongly 
nonlinear when the amplitude reaches a few tenths of mm, or more. At the same time 
the obtained coefficients of the series converge rapidly {l3/a ~ 0.15 mm^^), so that for 
1 2; I ^2 mm we can limit ourselves to expansion (|l]). Thus, Eq. ([T^ ) provides convenient 
analytical expression for the electric field in rather wide vicinity of the equilibrium particle 
position. 

Other experiments performed with different particle sizes and discharge conditions (Ar 
and Kr, p ~ 1 — 7 Pa) confirm that the sheath potential contains considerable anharmonic 
terms at low pressures. We have also investigated, how strongly the sheath potential is 
perturbed by the wire below the particle. Experiments with a wire at the same level with 





(14) 
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a particle (horizontal separation ~ 11 mm) show results, which are very similar to those 
reported here. 

In this Letter, we have presented one example of nonlinear particle oscillations in the 
sheath of a low-pressure rf discharge and its quantitative analysis. The proposed simple 
method opens a possibility to retrieve the profile of the electric field in the whole pre-sheath 
and upper part of the sheath region, using particles of different masses (and thus levitating 
at different heights). Analysis for each particle allows us to reconstruct the electric field 
in a range of a few mm, and therefore it is sufficient to use just a few different particles 
to get rather precise field distribution in the whole sheath range. Hence, we believe that 
experiments with excitation of nonlinear oscillations might be an effective way to study the 
sheath distribution. 
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FIGURES 



FIG. 1. Experimental setup. 

FIG. 2. Variation of the amplitude of particle oscillations close to the primary resonance for 
increasing (a) and decreasing (b) frequency of excitation, and for different magnitudes of the 
excitation voltage Uex- 50 mV (open circles), 100 mV (closed circles), and 200 mV (squares). Solid 
lines show the least-squares fit of the points using Eq. (P). The vertical dotted line indicates the 
position of the resonance frequency, ujq, obtained from the fit. 

FIG. 3. Variation of the amplitude of particle oscillations close to the superharmonic resonance 
for increasing (a) and decreasing (b) frequency of excitation, to, and for different magnitudes of the 
excitation voltage Uex- 4 V (open circles), 6 V (closed circles), and 7 V (squares). Solid lines show 
the least-squares fit of the points using Eq. (|To|) . The tilted dashed lines represent the amplitude 
offset [Eq. (^)]. The vertical dotted line indicates the position of the half resonance frequency, 
^ujQ, obtained from the fit. 

FIG. 4. Variation of the amplitude of particle oscillations close to the subharmonic resonance 
for increasing (a) and decreasing (b) frequency of excitation, uj. The excitation voltage C/ex = 2 V. 
Solid lines show the least-squares fit of the points using Eq. ([l3|) . The tilted dashed lines represent 
the amplitude offset [Eq. (0)]. The vertical dotted line indicates the position of the double 
resonance frequency, 2ujq, obtained from the fit. 
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